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Abstract 

We review recent progress towards the solution of exactly solved isotropic vertex 
models with arbitrary toroidal boundary conditions. Quantum space transformations 
make it possible the diagonalization of the corresponding transfer matrices by means 
of the quantum inverse scattering method. Explicit expressions for the eigenvalues and 
Bethe ansatz equations of the twisted isotropic spin chains based on the B n , D n and C n 
Lie algebras are presented. The applicability of this approach to the eight vertex model 
with non-diagonal twists is also discussed. 
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A vertex model is a classical statistical system defined on a lattice whose geometry is given 
by a possibly infinite set of straight lines on the plane pQ. The intersections of these lines are 
called vertices or sites and here we consider the situation where not more than two lines meet 
at every vertex. Clearly, the square lattice of size L x L is the simplest one and from now on 
we shall restrict ourselves to it. A physical state is then defined by assigning to each lattice 
edge a discrete variable having one out of q possible values. 

Next, we suppose that the corresponding row-to-row transfer matrix can be constructed 
from elementary local i-th site Boltzmann weights £43 (A) where A denotes a spectral parameter. 
This operator is best viewed as a q x q matrix on the auxiliary space A = C u whose elements 

L 

are operators acting on the Y[ <g>C? Hilbert space. Considering toroidal boundary conditions 

on the square lattice, the transfer matrix T(A) can be written in terms of the trace over A of 
the following ordered product of operator [3 01 El 

T(A) = Tr.4 [Q A C AL (X)C AL ^(X) . ..C A1 (X)} (1) 

where Q A are q x q c-number matrices, representing generalized periodic boundary conditions. 

A sufficient condition for integrability, i.e [T(X),T(fi)] = for arbitrary values of A and /i, 
is the existence of an invertible matrix -R(A, /1) satisfying the property [3 El 

R(X, n)C M (\) <g> C M (n) = C M {fi) ® C Ai (X)R(X, /1) (2) 

and that the matrix Q A is a possible representation, without spectral parameter dependence, 
of the quadratic algebra (J2J), namely |lj 

[R(X, ij),G a ® S a ] — 0, (3) 

As usual the i?-matrix -R(A, /x) is required to satisfy the famous Yang-Baxter equation 

-R23(Al, A2)-Rl2(Al, As)i?23(A2, A3) = Rl2{X2, A3)i?23(Al, A3)i?i2(Ai, A2), (4) 

In this paper we will consider integrable models whose corresponding i?-matrices are ad- 
ditive with respect the spectral parameters, R(X, fi) = R(X — //). In this case, the simplest 
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spectral parameter dependent representation of the Yang-Baxter algebra (0) turns out to be 

A* (A) = P M &(\), (5) 

where Pj^ is the exchange operator on the space A <S> Cf. 

If the matrix Qa is non-singular it is possible to derive an integrable quantum spin chain 
from T(A). Assuming that the operator is proportional to the permutator P^, say at 

certain special point A = 0, the corresponding one-dimensional Hamiltonian reads 



w _Vp ^Aa»(A) , fi _i p dC L1 (X) 

n ~ Z^^r - Ia=0 + ^ Pli ^a~ Ia=0 ^ l (6) 

When the boundary matrix is non-diagonal, the diagonalization of either the transfer 
matrix or the Hamiltonian (JHJ) is indeed a highly non-trivial problem in the field of integrable 
models. The main difficulty is concerned with the apparent lack of simple references states to 
start the Bethe ansatz analysis. Here we would like to present the steps towards the direction 
of solving integrable isotropic vertex models with non-diagonal toroidal boundary conditions. 
As concrete examples we will consider those systems whose rational i?-matrices are invariant 
by the B n , D n and C n symmetries. One way to construct rational solutions of the Yang-Baxter 
equation (£Q) is by means of the braid- monoid algebra [6 at its degenerated point [7]. This 
algebra is generated by the identity Ii, by a braid hi and a Temperley Lieb operator E{ acting on 
sites i of a chain of length L. On the degenerate point the braid operator becomes a generator 
of the symmetric group, namely 



z^e i} (7 ) 



a,b=l 

where are the q x q Weyl matrices acting on the space Cf. 

The monoid turns out to be represented by the following expression [H| 

a,b,c,d=l 

where the matrix a for the models B n , D n and C n are given by 



T -T I ( ^ nxn Znxn 

OtB n — - L 2n+lx2n+l, «D„ — - L 2nx2n, ®-C n I \?) 



-T O 



such that Tkxk is defined as a k x k anti-diagonal matrix. 

The set of algebraic relations satisfied by the braid and monoid at its generated point can 
be "Baxterized" in terms of rational functions. More specifically, the solution -R(A) in terms 
of combinations of the identity, braid and monoid is given by [7J |S] 

Ri, i+ i(\)=I i + Xb i -j^E i (10) 

where the values of parameter 5 are 

$B n = -n + l; 5 Cn = -n-l 5 Dn = -n + l (11) 

Let us now turn our attention to the diagonalization of the transfer matrix for the 
above B n , D n and C n vertex models, 1 considering the most general admissible boundary 
matrix satisfying the condition (|3*|) . i.e [Ei,GA® Ga} = 0. Denoting by Ma the matrix that 
diagonalize the boundary matrix Ga an d by inserting the terms M A Mj^ all over the trace (P) 
one we can write T(A) as 

T(A) = Tt a [M a D a M^M a (M^C al (\)M a )...(M a 1 Cai(\)Ma)M a 1 }, (12) 



Tr 



A 



DaCal{\)Cal-i{\)...Cax{\) , (13) 



where D A is diagonal matrix whose entries are the eigenvalues of Ga an d the ^-operators are 
given by 

C Al (\)=M A 1 C Al (\)M A . (14) 

The next step in our approach is to note that it is always possible to choose an invertible 
transformation Ui on the space Cf such that 

U^CAiiWi = C Ai (X) (15) 

and therefore to undo the modifications on the Lax operators ()14j) by means of quantum space 
transformations. 



1 We recall that similar problem for the A n model has been tackled in ref. 9 . 
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Now one can take advantage of this remarkable property by defining a new transfer matrix 
T'(A) 

L L 

T\X) = J] ®UJ 1 T{X) J] ®Uj = Tr A \D A C AL {X) . . . C A1 (X)}, (16) 
j=i j=i 

which is precisely the transfer matrix of the vertex model we have started with diagonal twist 
D A . 

Because the boundary D A is diagonal the transfer matrix T'(A) can be diagonalized with 
very little difference from the standard periodic case [Hj. Furthermore, the operators T(X) 
and T'(A) share the same eigenvalues and if is an eigenstate of T'(X) the corresponding 

L 

eigenvector of T(A) is then JJ®^!^'}- Considering that the algebraic framework to 

3=1 

diagonalize T'(A) has already been described in ref. [S], there is no need to repeat it here, and 
in what follows we shall present only the final results concerning the Bethe ansatz equations 
and eigenvalues of the related Hamiltonian (jHJ. The expression for the latter can be written 
in a compact form in terms of the underlying Cartan matrix C a b and the normalized length r] a 
of the roots. To each a-ath root we associate a set of rapidities \j that satisfy the following 
non-linear coupled equations, 



^(a) _ 6a 



2r la 



X " + 



m b \(a) \(b) C ayb 



_ II J = l,...,m.; a = l,...,n (17) 

V J 2 Va J ya+1 b=l k=l, k^j A j A k + 2 Va 

where g a is the a-th eigenvalue of the matrix D A . 

Before proceeding it should be remarked that due to the constraint [E iy Q A ® Q A ] = not 
all the eigenvalues g a are independent. It turns out that only the first n ratios are indeed 
arbitrary. The eigenvalues E(L) of the Hamiltonian ((HJ) are parameterized by the variables 
by 



E(L) = — — 7-- + L (18) 

tr t^i ] 2 - 1/4 

We expect that these results extend to all isotropic integrable vertex models invariant 
by the discrete representations of Lie algebras as well as to superalgebras. Note that such 
systems possess a broader class of possible non-diagonal boundary matrices as compared with 



their trigonometric counterparts . This means that isotropic vertex models with the most 
general twisted boundary conditions are in fact genuine systems that deserve to be studied 
independently. They are also of potential physical interest, since suitable combinations between 
non-diagonal boundaries and rational /^-operators can described interesting solvable atom- 
fields models jlOj . 

A natural question to be asked is whether or not this approach can also be of utility for 
non-rational vertex models. A tantalizing problem would be the solution of the eight-vertex 
model in the presence of non-diagonal toroidal boundary conditions. The symmetrical eigth- 
vertex model pQ possesses four different Boltzmann weights a(A), 6(A), c(A) and d(X) whose 
local operators C^i(X) which are given by the following 2x2 matrix 

C^-f^'^^ W + ) (19 ) 

V c(X)at + d{\)a- b{\)af a - + a(\)a7o+ J 

and af are Pauli matrices acting on the z-th sites of an one-dimensional lattice of size L. 
This vertex model is known to be solvable in the manifold 

= a\\) + b\\)-c\\)-d\\) = q(A)6(A)-c(A)d(A) 

o(A)6(A)+c(A)d(A) a(A)6(A) + c(A)d(A) 1 ' 

where A and Y are arbitrary constants. A possible non-diagonal twist compatible with inte- 
grability is given by 

One can now follow the same steps discussed above. Though we could not undo completely 
the modifications carried out on the auxiliary space due to the manipulations (|13|) we find out 
that the following quantum space transformation 

4*(A) = Ur 1 ^^ (22) 

where the matrix JJ{ is given by 

U t =\ | | (23) 



This transformation leads us to an operator C A i{X) that preserves the eight vertex form 
i~Ai{X) = _ (24) 

\ c(A)a+ + d(A)<7- 6(A)«7 i + t7- + a(A)a; a+ J 
where the new Boltzmann weights a(A), 6(A), c(A) and d(X) are given by 

_ = a(X) + b(X)+c(X) + d(X) 



6(A) = - ^ ^ - ^ (26) 



2 




a(A) + 6(A) - 


c(A) - tZ(A) 


2 




a(A) -6(A) + 


c(A) - d(X) 


2 




a(A) -6(A) - 


c(X)+d(X) 



c(A) = v ' ' v ' ^ (27) 



d(A) = ^>-^>-^>^^> (28) 

whose invariants are A = 4- and f = ■£ 2 . 

As a consequence of that our remaining task now consists in diagonalizing the following 
transfer matrix, 

f (A) = Tr^ [V A C AL (X)C M -i(X) ■ ■ • Aai(A)] (29) 
where the diagonal boundary matrix is 

Va = ( o -l ) (30) 

By construction T(A) (|29j) and T(A) given by Eqs. (jlJ19J21)) share the same eigenvalues while 
the eigenvectors are related by the similarity transformation (|23|). Though this procedure brings 
some simplification in the eigenvalue problem, it is not enough to make the diagonalization 

of the transfer matrix T(A) amenable by Bethe ansatz analysis. This is because the operator 
2 This then rccmphasize why the isotropic limit A = 1 is special under both auxiliary and quantum trans- 
formations. 
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A/u(A) has no simple local pseudovacuum that annihilate one of its off-diagonal matrix elements 
for arbitrary values of the spectral parameter. The standard way of solving this problem is by 
means of the so-called Baxter's gauge transformations PJ 0] which unfortunately does not work 
here since the diagonal boundary T>^ is not an identity matrix. In other words, the problem of 
finding gauge transformations Mj(A) with the conditions that both the transformed operator 
M^\(A)£^(A)Mi(A) has a local vacuum independent of A and that does not spoil the diagonal 
property of T>a has eluded us so far. An advantage of this approach, however, is that we can 
easily identify the existence of at least one case in each the eigenvalue problem for T(A) (J2HJ) 
can be solved by standard algebraic Bethe ansatz. This clearly occurs when the Boltzmann 
weight d(X) (|28|) is null. Direct inspection reveals us that this happens at the point in which 
the modulus /C of the elliptic functions parameterizing the eight vertex Boltzmann weights 
becomes unity. More specifically, at the value /C = 1 the weights a(A), b(X), c(A) and d(X) are 
given by the following expressions 

a(A) = tanh[A + 7] b(X) = tanh[A] (31) 



c(A) = tanh[7] d(X) = tanhfy] tanh[A] tanh[A + 7] (32) 

which due to (}2"%j) implies d(X) = 0. 

Thanks to the above simplification it now remains only the diagonalization of a symmetric 
six vertex model, whose solution has appeared in many different contexts in the literature. 
The result for the eigenvalue A(A) of T(A) is therefore 

where the weights a(A), 6(A) and c(A) are 

3 mh[A + 7 ] sinh[Al sinhhj 

cosh [A J cosh [7 J cosh [A + 7J cosh [7 J cosh [A + 7J cosh[AJ 

The integers m < L parameterize the multiparticle state of T(A) and the Bethe ansatz 
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roots Xj satisfy the equations 



rn 



b(X k - Xj) a(\j - X k ) 
a(X k - Xj) b(Xj - X k ) 



n 



3 = 1, • • • , 



m 



(35) 



It is conceivable that an adaptation of the above ideas might work for the general eight 
vertex model (fTHj) with the boundary (j21j) . In fact, the case K, = 3 was solved by means 
of certain functional relations even though the eigenvectors structure is not yet known. If this 
could be carried out, even for particular values of the modulus /C, it would be an important 
step toward the understanding of properties of the eight vertex. 
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